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. Abstract 

■ A general formulation for describing odd-harmonic cosmic strings is 
I developed and used to determine the self-intersection properties of high- 

■ harmonic loops. This is important because loop formation mechanisms 
0^ I produce high-harmonic components (kinks) which can only be eliminated 

■ very slowly by gravitational radiation, damping by the dense surrounding 
C^hI plasma in the era of string formation, or by the expansion of the Universe. 
^ • For the class of loops examined it has been found that in the high-harmonic 
Qj ' limit, essentially all cosmic loops self-intersect. 
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1 Introduction 

Topological defects such as monopoles, cosmic strings and textures, are pre- 
dicted by a number of grand unified theories as a result of symmetry breaking 
and may have important comological consequences [p. Cosmic strings, in par- 
ticular, could be responsible for the creation of the density fluctuations needed 
for large-scale structure formation |2|, ^. They are characterized by a tension 
or mass per unit length fi which is related to the temperature Tqut at which 
the phase transition takes place by /x ~ ^gut- Cosmic strings can meet and 
intercommute (exchange partners) or self-intersect. The latter process results 
in the formation of cosmic loops. 

In the old cosmic string scenario |Q, |5|, ^ it was envisaged that these loops 
gave rise to the density fluctuations responsible for galaxy, galaxy cluster and 
supercluster formation. Loops decay primarily via gravitational radiation, giv- 
ing them a very long lifetime. Recent numerical simulations of cosmic string 
networks |Q, ^, 10 1, while confirming that a scaling solution is probably 



reached, have shown that there is a great deal of small-scale structure on the 
strings, and that the typical size of the loops formed is very small compared to 
the horizon. Since the lifetime of a loop is proportional to its length, this means 
that they do not live very long, and are unable to act as the seeds around which 
galaxies form. The dominant density perturbations must come from the long 
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strings. Nevertheless the loops could still have significant effects collectively. 
The extent to which they do so depends strongly on whether they typically 
break up into large numbers of small loops or reach a stable form that can 
survive a long time. 

It is difficult to analyze the fate of such loops using the numerical simulations 
of cosmic string networks, because they lack the necessary resolution. Bennett 
and Bouchet ^ for example estimated that most loops break up into about ten 
stable daughter loops, but since the resulting stable loops are comparable in 
size to the resolution, it is unclear how far we should trust this conclusion. So 
it is desirable to study the evolution of loops analytically. 

Several string parametrizations, corresponding to low harmonics (first and 
third), have been analyzed for self- intersections — those of Kibble and Turok 
a one-parameter family of cosmic loops, Turok with two parame- 
ters, Chen, DiCarlo and Hotes [|l3|, with three and finally DeLaney, Engle and 



Scheick |14|, with five. The principal justification for limiting study to such 
configurations (aside from computational convenience) is the argument that 
high-harmonic components are damped either by the expansion of the Uni- 
verse, by friction with the dense surrounding plasma in the era immediately 
after string formation, or by the gravitational back-reaction, so that after a 
time loops become relatively smooth [jl2|, |6|, 14]. However, we shall argue 
that these mechanisms may not in fact be so relevant to loop stability. 



Figure 1: String intercommutation (a) and loop formation mechanisms: By the 
self-intersection of a long string (b) and by the self-intersection of a cosmic loop 
(c). 

Loop formation mechanisms are illustrated in Figure 1. When strings 
meet they intercommute (Fig. 1(a)). Cosmic loops can be created by the self- 
intersection of a pre-existing loop (Fig. 1(b)) or of a long string (Fig. 1 (b) ) . Both 
mechanisms inevitably produce kinks, or equivalently, high-harmonic compo- 
nents. Damping by interaction with the dense surrounding plasma only affects 
strings shortly after their initial formation. However, the loops of cosmological 
interest are produced when these effects are no longer important. Similarly, 
the expansion of the Universe is relevant only for loops comparable in size to 
the horizon, not to the small loops that dominate the recent loop distribution. 
High harmonics can also be radiated away by the emission of gravitational 
waves. The time scale of gravitational radiation effects, however, is very much 
larger than the period of oscillation of a cosmic loop. Therefore, since a cosmic 
loop in a self-intersecting trajectory will self-intersect within one oscillation pe- 
riod, it is reasonable to consider gravitational effects negligible as far as loop 
stability is concerned. 

For all these reasons, it would seem therefore to be useful to investigate 
the properties of higher-harmonic cosmic loops to determine how the self- 
intersection probability varies with the number of harmonics. This is the prob- 
lem we address in this work. The presence of these extra harmonics makes 
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the parameter space of cosmic loops larger and more complex than that of the 
string configurations analyzed so far. 

Section 2 provides a summary of the string equations of motion and the gen- 
eral solutions in terms of Fourier series and in terms of the product of rotations 
1 15]. Section 3 develops a general formulation for an odd-harmonic cosmic string 
in terms of the product representation, and describes our method for finding 
self-intersections. An explanation of how it is computationally implemented is 
given in an Appendix. In section 4 we describe the self-intersection results ob- 
tained for harmonic parametrizations which range from a 3/3 string {i.e., one 
including up to the third harmonic on each of the left- and right-moving halves) 
up to a 41/41 string. We also discuss their sensitivity to different cutoffs. Fi- 
nally, Section 5 contains a discussion of the new self-intersection results, their 
effect on galaxy-formation scenarios, and new directions of research. 



2 Equations of Motion and General Solutions 

The position of a string is described by x^{cr, r), where a and r are the spacelike 
and timelike variables parametrizing the world sheet that the string sweeps out 
in space-time (c is taken to be 1). It is convenient to use the orthonormal (or 
conformal) gauge defined by the constraints 1 1£ ] 

daX^drXfj, = daX^daX^ + drX'^ drX ^ = 0. (1) 

For loops small compared with the horizon size and when gravitational radia- 
tion and external gravitational fields are considered negligible, the equation of 
motion of a cosmic string is the classical relativistic wave equation 

dlxi" - a^x'^ = 0. (2) 

We are also free to impose the further condition x° = t, so that x^{a,t) = 
{t,r{a,t)). Then the solution in the centre of mass frame can be written 



ri(T,t) = -[^iu)+hiv)], 



(3) 
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These vectors trace out closed loops on a unit sphere [|lT|. Here they will 
be referred to as the two halves of the string. 

We can choose, for convenience, the period in a to be 27r. Thus r is periodic 
in t with effective period it. The two halves of a closed string loop can be 
expanded in terms of the Fourier series 



a' = Z + X) A„ cos (nn) + X] B„ sin(nti), 
b' = — Z + C„ cos(nt') + D„ sin(nt') 



(5) 
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where a' and b' do not, in general, have the same number of harmonics. The 
hnear centre of mass (cm.) terms of a and b must be related due to spatial 
periodicity of closed loops. The problem therefore is reduced to finding a set of 
harmonic coefficients that satisfy the constraint equations (^) and the relation 
between the cm. terms of a and b (^). 

The general solution of an A^-harmonic string half requires the satisfaction 
of 4:N + 1 nonlinear relations between the vector coefficients. For the A^=l, 
N=2 and A^=3 cases these relations can be solved with relative ease but for 
higher harmonics they become intractable. An alternative procedure is available 
involving products of rotations |15|. In terms of this representation a general 
A^-harmonic string half is given by 



c^('u) = pn+iRz{u)pnRz{u) . ■ ■ p3Rz{u)p2Rziu)pii. 



where 



Pi = p{o^,(i)i) = Rz{-ei)R,{(t)i)Rz{ei 

and the rotation matrices are 



(6) 



(7) 
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and 



RM 



cos I 
sin ( 



— sm ( 
cos 6 



It is clear that the magnitude of the vector remains fixed because all that is ever 
done is to rotate the original vector a certain number of times. The product 
representation is also complete but the proof of this is more involved and will 



not be included here for brevity (see Ref. |15l). It can be seen, however, that it 
exhibits the correct number of degrees of freedom. A general A^-harmonic string 
half in terms of the Fourier series has 6A^ + 3 vector coefficients (Eq. ^) which 
must satisfy 4A^ + 1 nonlinear relations (Eq. leaving a total of 2N + 2 inde- 
pendent degrees of freedom (Eq. ^) . The ranges of the independent parameters 
9 and (j) are polar-like and azimuthal-like respectively. 



3 The General Odd-Harmonic String 

The magnitude constraint has been solved by generating the string using the 
product of rotations. The relation between the centre of mass terms in the two 
string halves due to overall spatial periodicity of closed loops can be solved by 
taking an odd-harmonic string with no zeroth harmonic component. It should 
be noted that there are many other harmonic parametrizations with no zeroth 
harmonic. We have chosen the odd-harmonic string because it is the most 
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simple configuration of this type. With this procedure there will be a loss of 
generality — the loops have an extra inversion symmetry — but the problem 
of finding a set of parameters which satisfies the relation between the centre 
of mass terms of any two string halves in terms of the product representation 
resists analytical solution. 

The odd A^-harmonic string half is given by 

c'Niu) = pN+2Rzi'2u)pNRz{'^u) . . . Rz{2u) psR^iu) piZ, (8) 

where N is an odd integer. Because we want no zeroth harmonic component 
in the string we must take the angle 0i = it^ in both halves of the string. 
This gives a total of four families of strings per harmonic configuration just 



as found in Ref. |14] for the general three harmonic string. In this work we 
will only analyze those with = ^ in both halves of the string. For an N/M 
odd-harmonic string (the first letter denotes the number of harmonics of a and 
the second that of b) this leaves us with a total of iV + M + 4 free parameters 
((A^-|-2)-|-(M-|-2)). Two more parameters, the 9i for a' and b', control the origin 
of u and v respectively and are set to — ^ to put the vector into standard form 



1 15]. Because we are only looking for self-intersections we can eliminate three 
extra parameters from overall orientation freedom, namely 9^-^-2 from a' and 
9n+2 and (pN+2 from b' leaving Rx{4)n+2) in a' to ensure relative orientation 
freedom. This leaves a grand total of M -|- — 1 free parameters for the most 
general N/M odd harmonic cosmic loop with no zeroth harmonic. This result 
agrees with the number of string parameters of the 3/1 string presented in |13| 
and the 3/3 string presented in 1 14 1 . 

When all the redundant parameters have been eliminated the harmonic 
expansions in terms of the product of rotations for the two string halves are 

a^(n) = Rx{(I)n+2)Rz{2u)pnRz{2u) . . . R,{2u)pMu)^ (9) 

and 

h\j{v) = R,{2v)pmRz{2v) . . . R,{2v)pMv)5c (10) 

These formulae are the expressions for the derivatives of the string halves. 
What we need in order to find self-intersections, however, are the coefficients 
of the string halves themselves. Fortunately a recurrence relation developed in 
Ref. ||l5|| can be used to compute the vector coefficients of the Fourier series in 
terms of the angles 9i and It yields the vector coefficients in standard form, 
so the coefficients of a still need to be multiplied by Rxi4'N+2)- 

We have thus found a means of generating the vector coefficients of a general 
N/M odd-harmonic string which satisfies the magnitude and centre of mass 
constraints. This is the set of strings which we will analyze for self-intersections. 



(a) (b) 



Figure 2: Two examples of strings generated by the product of rotations and 
the recurrence relation: (a) A 19-19 harmonic string and (b) a 47-47 harmonic 
string. 
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The equation to be solved for finding self-intersections is 



r{a,t) = ria',t) (11) 
wliich in terms of the Fourier series is 

— [A„ sin(nn) — B„ cos(nn) + Cn sin(nf) — D,^ cos(nf)] 

= ^ — [A„ sin(nii') — B„ cos(nti') + Cn sin(n?j') — D„ cos(nti')] (12) 

where 

u = a — t, V = (7 + t, (13) 

u = a' — t, v' = a' + 1, (14) 

and a a' + 2mr for any integer n. This yields three non-linear equations in 
a, a' and t to be solved simultaneously. Although individual cases of string 
parametrizations [11, |l^, |l^, 17] have been solved analytically this problem 



is intractable for a general N/M odd-harmonic string. Numerically, however, it 
can be solved. The general procedure used here is described in the Appendix. 
To eliminate spurious intersections, it is necessary to introduce a small-distance 
cutoff. Self-intersection results are typically [^, 0] given as a function of the 
cutoff, and comparisons between them should only be made for the same cutoffs. 



4 Self-Intersection Results 

The Chen, DiCarlo, Hotes (CDH) string |l^ has been used to test the cor- 
rectness of the self-intersection method developed here. The percentage of self- 
intersections for the CDH string has been found to be 33.6% it 5%. This result 
is in very good agreement with all the previous work ||l3| , p^ . The Turok 
string [|l^, 13] has also been tested for self- intersections and a percentage of 



2% lb 1.25% was found, which is in excellent agreement with the previous work 
as well. Furthermore, the 3/3 string presented here, which is equivalent to the 
DeLaney, Engle and Scheick (DES) string [0] has an intersection probability of 
about 0.6 (see Figures 3, 4 and 5), this is also consistent with their results. The 
self-intersection probabilities obtained for the generalized odd-harmonic strings 
are presented in the figures below. The parametrizations for which intersec- 
tions have been calculated range from the 3/3 string to the 25/25 string for the 
smallest and largest cutoffs and from the 3/3 string to the 41/41 string for the 
middle cutoff. 

A random flat parameter distribution has been used for the angles in the 
vector coefficients and the number of points K along the string was chosen to 
he K = 600, yielding a resolution r] ~ 0.0104712 radians which appears to 
be sufficient. Three different cutoffs have been used, 6a = \a — a'\ = 0.084 
radians, 0.126 radians and 0.168 radians, corresponding to eight, twelve and 
sixteen step lengths, respectively. These are approximately the lower, middle 
and upper cutoffs used in the previous work on self-intersections thus allowing 
easy comparison. Each data point on the graphs is an average value of ten 
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Figure 3: Self-intersection probability as a function of the number of harmonics 
for a cutoff of Sa = \a — a'\ = 0.084 radians. 



Figure 4: Self-intersection probability as a function of the number of harmonics 
for a cutoff of 6a = \a — a'\ = 0.126 radians. 

samples of a hundred strings each. The errors are the standard deviation from 
the average of these ten samples. 

In Figures 3, 4, and 5 plots are shown of the self-intersection probability as 
a function of the number of harmonics for the three cutoffs. It should be noted 
that the harmonic parametrization corresponding to a value N on the x-axis on 
the figures corresponds to a N/N harmonic string. The only parametrizations 
analyzed in this work correspond to strings with the same number of harmonics 
in each half. 

In Figure 6 a plot of the logarithm of the stability probability as a function 
of the number of harmonics for the three cutoffs is presented. The stability 
probability is defined as the probability that loop will not self intersect (i.e. 
Pstab = 1 — PSi)- The last six points do not have errorbars due to the fact that 
the lower value of the error is at infinity. As can be seen the stability probability 
Pstab can be fitted quite well to an exponential curve. The dashed curve on the 
plot is given by 

Pstab = e"+^^ (15) 
with parameters a = — 0.4 and (3 = —0.2. The continuous curve is 

Pstab = iV^e^+'^^ (16) 

with parameters S = —0.2, 7 = —0.4 and k = —0.14, where is the number of 
harmonics. It is clear that in the high harmonic limit all loops will self-intersect. 



5 Discussion and Conclusions 

We have set out to examine the self-intersection properties of high-harmonic 
cosmic loops. After the very early stages of string formation, damping by the 
dense surrounding plasma, by the expansion of the Universe, or by the emission 
of gravitational radiation does not have significant effects on the time scale of 



Figure 5: Self- intersection probability as a function of the number of harmonics 
for a cutoff oi da = \a — a'\ = 0.168 radians. 
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Figure 6: Plot of the logarithm of the stability probability as a function of the 
number of harmonics for the three cutoffs. The strings plotted here are the 
same as in the previous figures. 



the order of the oscillation period. Consequently, these processes are irrelevant 
to the question of how many loops self-intersect. 

We have found the most general N/M odd-harmonic string in terms of the 
product of rotations and devised a means for finding self-intersections. The 
results agreed with all previous work on the subject, but extended it to much 
higher harmonics. We have systematically analyzed the resulting string con- 
figurations to find the intersection probability in the high-harmonic limit. We 
have found that as a function of the number of harmonics the intersection prob- 
ability is well described by an exponential. In the high-harmonic limit almost 
all loops self- intersect. 

Although loops have been found not to be as important for structure for- 
mation in recent numerical simulations of cosmic string networks as was once 
thought, they may still have significant effects, on the density perturbations Q, 
on the microwave background [18|, on the gravity- wave background and 
possibly on the baryon asymmetry |20|. On the basis of our results, it appears 
that cosmic loops formed after the initial string formation era by the mecha- 
nisms described in Section 1 will almost all self-intersect within the first period 
of oscillation. As a result of the intersection process two or more kinky daughter 
loops are produced and, by the same line of reasoning, it is extremely probable 
that the daughter loops will re-intersect. It is not certain how far this process 
will go, but it seems likely to continue until the resulting loops have only a very 
few kinks, by which time they will be extremely small. Clearly it would be 
desirable to follow the evolution of the loops through several generations to see 
whether they continue to divide. Work in this direction is in progress. Bennett 
and Bouchet ^ and Allen and Shellard [0] found in their simulations that 
the size of most of their child loops was determined by the resolution of their 
simulations, i.e., close enough to the lower cutoff on loop size to suppress the 
chances of further division. In this work we have formulated loop trajectories 
analytically in an attempt to complement the numerical simulations made so 
far. We believe our method can provide a more definite answer to the question 
of the fate of cosmic string loops. 

At birth, the loops are already fairly small compared to the horizon. A 
sequence of self-intersection processes will turn the initial parent loop into tiny 
loops, with very short lifetimes, so that they will decay into relativistic particles 
within a relatively short period of time. This is likely to reduce many of the 
observable effects of loops, in particular their effects on density perturbations, 
on the microwave background and on the gravitational-wave background. But 
it might enhance the effect on the baryon asymmetry. 
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Appendix 



Here we desribe the implementation of our analysis. 

When the number of harmonics for each half of the string has been selected, 
a set of string coefficients is randomly generated using the recurrence relation 
mentioned in Section 3. The number K of points along the string must be 
chosen. It determines the resolution r] of the string in a and t, via r/ = 2tt/K. 
For consistency the resolution in a and t should be the same. It should be small 
enough to ensure that the small scale detail of the loop is not lost. 

The string is then time evolved for the given number of timesteps K/2, with 
t ranging from to vr, the effective period. At each timestep ti, we search for 
self-intersections by expanding the equations (^2|) to first order about every pair 
of points Uj, Vj and n^, given by 



<7fc — ti, 



Gj + ti, 
+ti. 



where j ranges]^ from 1 to K/2, and k from j to K/2 and from j + K/2 to K, 
as follows 

— I A„[sin(nnj) + ncos{nuj){u — Uj)\ — Bn,[cos(nnj) — nsm.{nuj){u — Uj)] 

+C„[sin(ni;j) + ncos{nvj){v — Vj)] — D„[cos(nfj) — nsui{nvj){v — fj)]} 

= — |A„[sin(ntifc) +71 cos(ntifc)(u' — Ufc)] — B„ [cos (nn^) — nsva.{nuk){u' — Uk)\ 

+C„[sin(nt;fc) + ncos{nvk){v' - Vk)] - T>n[cos{nvk) - nsm{nvk){v' - 
Substituting for u, v, u' , and v' from ([l4[ ) and rearranging the terms yields 

o" ^1 A„ cos(ntij) + B„ sin(nnj) + C„ cos{nvj) + D„ sin(nf j)| 

-o-'^|a„ cos(nufc) +B„ sin(nnA;) + C„cos(ntife) + D„ sin(nwfc)| 
+t A„[cos(nn-,) — cos{nuk)] — Bn[sin(nnj) — sin(nnA;)] 

+C,„[cos(nt'j) — cos{nvk)] + D„[sin(nfj) — sin(nt'fc)]| 



— sm(nuj) — Ui cos(nuj) 



n 



— smlnuk J 
n 



Uj cos[nuk, 



-Br 



+ Cn 



n 



cos{nuj) + Uj sm{nuj] 



1 



— sin(nt;j) — Vj cos{nvj) 



n 



cos{nuk) + Uk sm(nuk) 



— sm{nvk) — Vk cos{nvk) 



^The symmetry of this set of strings is such that an intersection found at a point ai on 
the string will mean that there is another intersection at (Ti + vr, so we only need to search for 
self-intersections in half the total range in a. 
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^— cos{nvj) + Vi sm{nvj)^ — ^— cos{nvk) + Vj sm(nvk] 



As can be seen this is a linear system of three equations and three unknowns 
a, a' and t, which is readily solvable by Cramer's rule. An intersection is found 
when the conditions \ti—t\ < r//2, C7| < ry/2 and \ak—a'\ < r//2 are satisfied. 
This makes the choice of resolution crucial to the accuracy of our results. The 
process is systematically repeated for every timestep until an intersection is 
found. 

To avoid trivial solutions (i.e., a = a') a cutoff Sa has been introduced. 

This means that if an intersection occurs at aj and ak with \aj — ak\ smaller 
than the cutoff it is neglected. In practice, the Taylor expansion is not made 
for pairs of points for which \aj — (Tfc| < 5a. 
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